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ABSTRACT 

Large samples of high-redshift supernovae (SNe) are potentially powerful probes 
of cosmic star formation, metal enrichment, and SN physics. We present initial results 
from a new deep SN survey, based on re-imaging in the R, i', z' bands, of the 0.25 deg 2 
Subaru Deep Field (SDF), with the 8.2 m Subaru telescope and Suprime-Cam. In 
a single new epoch consisting of two nights of observations, we have discovered 33 
candidate SNe, down to a z'-band magnitude of 26.3 (AB). We have measured the 
photometric redshifts of the SN host galaxies, obtained Keck spectroscopic redshifts 
for 17 of the host galaxies, and classified the SNe using the Bayesian photometric 
algorithm of Poznanski et al. (2007) that relies on template matching. After correcting 
for biases in the classification, 55% of our sample consists of Type la supernovae and 
45% of core-collapse SNe. The redshift distribution of the SNe la reaches z ~ 1.6, 
with a median of z « 1.2. The core-collapse SNe reach z w 1.0, with a median of 
z w 0.5. Our SN sample is comparable to the Hubble Space Telescope / 'GOODS sample 
both in size and redshift range. The redshift distributions of the SNe in the SDF and 
in GOODS arc consistent, but there is a trend (which requires confirmation using 
a larger sample) for more high-z SNe la in the SDF. This trend is also apparent 
when comparing the SN la rates we derive to those based on GOODS data. Our 
results suggest a fairly constant rate at high redshift that could be tracking the star- 
formation rate. Additional epochs on this field, already being obtained, will enlarge 
our SN sample to the hundreds, and determine whether or not there is a decline in 
the SN la rate at z > 1. 

Key words: supernovae: general — cosmology: observations, miscellaneous — sur- 
veys 



1 INTRODUCTION 
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Supernovae (SNe) hold the key to several open astrophysical 
questions. In particular, SNe are the main source of heavy 
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elements, and their energy input to the interstellar medium 
(ISM) is a vital ingredient in galaxy formation (e.g., Yepes 
et al. 1997). Learning the nature of the progenitors of dif- 
ferent SN types, the SN rates, and the SN distributions in 
both space and cosmic time, are essential steps toward un- 
derstanding metal enrichment and galaxy formation (e.g., 
Kobayashi, Tsujimoto & Nomoto 2000), and will be a prime 
topic of study with the James Webb Space Telescope. 

It is widely believed that SNe la arise from the explo- 
sions of white dwarfs (WDs) in binary systems, while all 
other SN types result from massive-star core collapse (e.g., 
Filippenko 1997; Gal- Yam et al. 2007; Li et al. 2007, and 
references therein). Core-collapse SNe (CC-SNe) therefore 
explode promptly (< 10 Myr) after the formation of a stellar 
population, and the CC-SN rate will trace the star-formation 
history (SFH). In contrast, SNe la should occur only after 
WD formation and binary evolution, with a "delay" of order 
0.1-10 Gyr. The SN la rate vs. cosmic time will therefore 
be a convolution of the SFH with a "delay function" — the 
SN la rate following a brief burst of star formation. 

The progenitors of SNe la have not been identified, and 
could be WDs accreting from main-sequence companions, 
sub-giants, or giants, or alternatively "double-degenerate" 
WD pairs that merge (e.g., Hillebrandt & Niemeyer 2000, 
and references therein). Each of these scenarios predicts a 
different delay function (e.g., Kobayashi et al. 1998; Madau, 
Delia Valle & Panagia 1998; Yungelson & Livio 1998; Bel- 
czynski, Bulik & Ruiter 2005; Greggio 2005; see recent up- 
dates in Panagia, Delia Valle & Mannucci 2007), which in 
turn dictates the SN rate vs. redshift. For the local Universe, 
current estimates of SN rates (e.g., Cappellaro, Evans & Tu- 
ratto 1999; Mannucci et al. 2005; Sharon et al. 2007) are still 
uncertain, but improving substantially (Leaman, Li, & Fil- 
ippenko 2007, in preparation). At z w 0.3— 1.6, a number of 
sometimes conflicting measurements exist (Pain et al. 2002; 
Gal- Yam, Maoz & Sharon 2002; Tonry et al. 2003; Dahlen 
et al. 2004; Cappellaro et al. 2005; Barris & Tonry 2006; 
Neill et al. 2006; Sullivan et al. 2006b). 

Recent attempts to constrain the cosmic SFH and the 
SN la delay function using SN rates have yielded some in- 
triguing results. Gal- Yam et al. (2002) carried out the first 
measurement of cluster SN la rates at z = 0.25 and z = 0.9 
using deep (R w 27 mag) archival Hubble Space Telescope 
(HST) images of eight massive clusters. In Maoz & Gal- 
Yam (2004), the measured rates were compared to predic- 
tions normalized by the requirement that SNe la produce 
the observed mass of iron in clusters. The low observed rate 
atzwl implies that the SNe needed to produce the iron 
exploded at earlier times, arguing against models with long 
SN la time delays, which predict SN rates ten times higher 
than observed. Thus, cluster iron production by SNe la ap- 
peared to be a viable option only if SNe la have relatively 
short (< 2 Gyr) time delays. In contrast, comparison of the 
SN la rate vs. z in the field to the cosmic SFH has generally 
suggested a long delay time (Pain et al. 2002; Tonry et al. 
2003; Gal- Yam & Maoz 2004; Strolger et al. 2004; see, how- 
ever, Barris & Tonry 2006, who find evidence for a short 
delay). These results may therefore point to CC-SNe from 
an early generation of star formation with a top-heavy initial 
mass function, rather than SNe la, as the source of cluster 
enrichment. 

However, a re-analysis by Forster et al. (2006) of 



the iTST-based Great Observatories Origins Deep Survey 
(GOODS) data of Strolger et al. (2004), originally obtained 
to measure the expected early deceleration of the Universe 
(Riess et al. 2004a), indicates that these conclusions may be 
strongly dependent on the assumed SFH. Furthermore, evi- 
dence has surfaced that there may be two separate SN la 
channels, a "prompt" (~ 1 Gyr) channel and a "tardy" 
(~ 10 Gyr) one (Mannucci et al. 2005; Scannapieco & Bild- 
sten 2005; Mannucci, Delia Valle & Panagia 2006). A single 
delay distribution that is peaked at short delays but includes 
a long tail extending to large delays may be able to repro- 
duce most of the current data. Nevertheless, some rate mea- 
surements at similar redshifts are discrepant by large factors 
(e.g., Pain et al. 2002 vs. Barris & Tonry 2006). There is on- 
going debate as to whether incompleteness or contamination 
of SN samples are affecting the low or high-rate measure- 
ments, respectively. The decrease in the field SN la rate at 
z > 1 measured by Dahlen et al. (2004) and Strolger et al. 
(2004), if real, does not permit the existence of the prompt 
SN la component indicated by the high SN la rate seen in 
star-forming galaxies (Mannucci et al. 2005, 2006; Sullivan 
ct al. 2006b). Scannapieco & Bildsten (2005) have actually 
predicted that this apparent decrease will disappear with 
improved data. Barris & Tonry (2006) have argued that the 
claims for long delay times depend strongly on the pres- 
ence of very few SNe at the highest redshifts in the GOODS 
survey. Indeed, small-number statistics are a major limiting 
factor in all current estimates of SN rates at high z. Progress 
in this field requires larger high- z SN samples than presently 
available. 

However, with typical magnitudes of ~25-26, samples 
of high-z SNe are currently difficult to confirm (or classify) 
spectroscopically. Photometric classification of SNe has been 
shown do be a viable option by Poznanski et al. (2002), 
and further explored since then (e.g., Gal- Yam et al. 2004b; 
Riess et al. 2004b; Sullivan et al. 2006a). Recently, Poznan- 
ski, Maoz & Gal- Yam (2007, hereafter P07) have developed 
an automatic Bayesian algorithm for the purpose of classify- 
ing SNe with single-epoch photometry. The major novelties 
in the approach of P07 lie in the probabilistic classification, 
and the ability to quantify the success fractions, uncertain- 
ties, and biases of photometric SN typing. This opens the 
possibility for surveys for large numbers of SNe, and their 
classification, using solely single-epoch, ground-based, imag- 
ing. 

In this paper, we present initial results from a survey 
in which we attempt this new approach for obtaining large 
high-z SN samples. The survey is conducted by re-imaging 
the Subaru Deep Field (SDF; Kashikawa et al. 2004, here- 
after K04) in three bands. Keck spectroscopy of over half 
of the SN hosts, and for several hundred galaxies similar to 
the hosts, is used to estimate the reliability of our photo- 
metric redshifts. The photometry of the SNe and their hosts 
is used to classify the SNe with the algorithm of P07. As 
we show below, we identify 33 SNe, classified as either Type 
la or core collapse, with redshifts out to ~ 1.7, a number 
and redshift range comparable to those in the ffi>T/GOODS 
sample. While the redshift distributions of the samples are 
consistent, there is a trend for more high-z SNe la in the 
SDF. We derive the SN la rate, in four redshift bins, and 
find that while it is consistent, to within errors, with previ- 
ous results, at redshift z « 1.6 it is about 50% higher than 
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the value determined based on GOODS data, and at redshift 
z ~ 0.7 the rate we find is lower than that of GOODS. We 
are currently in the process of increasing our SN sample by 
obtaining additional epochs of the SDF, with the objective 
of acquiring a sample of several hundred high-z SNe. Those 
observations and their analysis will be described in future 
papers. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Imaging 

Our SN survey is carried out in the SDF (a = 13 h 24 m 39 s , 
S = +27°29'26"; J2000). The reference images for this field, 
obtained with the Suprime-Cam camera on the Subaru 8.2- 
m telescope on Mauna Kea, Hawaii, are described by K04. 
Suprime-Cam is a 5 x 2 mosaic of 2k x 4k CCDs at the 
prime focus of the telescope, with a field of view of 34' x 27', 
and a scale of 0'.'202 pixel -1 . This reference epoch of the 
SDF was obtained by deep imaging, in five broad-band filters 
(B,V,R,i', and z') and two narrow-band filters (JVB816 
and A r B921), of an area of 30' x 37', down to 3<r limiting 
magnitudes of B = 28.45, V = 27.74, R = 27.80, i' = 
27.43, and z = 26.62 (5a limits of B = 28.06, V = 27.33, 
R = 27.37, i' = 27.13, and z = 26.32; here and throughout 
the paper, all magnitudes are on the AB system; Oke & 
Gunn 1983). The imaging for the reference SDF epoch was 
carried out between April 2002 and April 2003, with some 
preliminary data collected in 2001, as described by K04. 

We re-imaged the field on 2005 March 5 and 6 (UT 
dates are used throughout this paper), in the three reddest 
Suprime-Cam broad bands: R, i , and z . These bands were 
selected as being the most efficient for finding high-redshift 
SNe, and for later classification (e.g., Poznanski et al. 2002; 
Gal- Yam et al. 2004b; Riess et al. 2004b). We obtained 22 
exposures of 360 s each in R, 36 of 300 s each in i' , and 76 
of 240 s each in z' , for a total of 7920 s (R), 10,800 s (i'), 
and 18,240 s (z'). We followed a dithering pattern similar 
to the one described by K04. Our seeing throughout the 
observations averaged 0'.'65 in all three bands, ranging from 
0'.'5 to 0'.'8. 

We also observed the field with the MOSAIC-I wide- 
field imager on the Kitt Peak National Observatory (KPNO) 
4-m Mayall telescope, on 2005 March 2, April 5, and June 
8, and on 2006 June 19. The MOSAIC-I camera consists of 
eight CCDs with a pixel size of 0'.'258 (Muller et al. 1998), 
covering a field of view of 36' x 36'. Each observation was 
in white light (i.e., unfiltered), typically 1.5-2 hr long, split 
into dithered subexposures, with a median seeing of about 
0'.'8. Compared to the R band, the unfiltered observations in- 
crease the throughput of stellar sources, depending on their 
colours, by a factor of 2-4, and also increase the level of 
the sky background by a factor of 4-8, depending on con- 
ditions. The signal-to-noise ratio (S/N) is therefore multi- 
plied by (2 - 4) / V4 - 8 w (0.7 - 2). Consequently, in many 
cases, unfiltered observations permit obtaining a deeper im- 
age, closer to the depth reached with the larger-aperture 
Subaru telescope, albeit with no colour information. The 4- 
m data provide us with monitoring over longer times, which 
is important for identifying variable active galactic nuclei 
(AGNs), the main contaminant in our SN search (see fj3|. 



The HST archive includes Advanced Camera for Sur- 
veys (ACS) images of a small portion of the SDF, which we 
use to examine the morphology of one SN host galaxy that 
lies in this field. The data were obtained with HST/ ACS in 
2002, on five separate occasions between May 20 and July 
3, in the F850LP band as part of program GO-9075 (PI: S. 
Perlmutter). The total exposure time is 5370 s. 

Subaru data were reduced following K04, using the 
Suprime-Cam pipeline SDFRED (Yagi et al. 2002; Ouchi 
et al. 2004), and involving the following main steps. The in- 
dividual frames were overscan subtracted, flat fielded using 
superflats, distortion corrected, matched to a common point- 
spread function (PSF) of 0'.'98, sky subtracted, registered, 
and combined. The combined image was then matched to 
the reference image by using the astrometri^ code to find 
the astrometric correction, and the IRAF (Tody 1986) task 
wregister was used to register the two images. Inducing PSF 
degradation to the new images may appear questionable, 
since it reduces the frame depth and could adversely af- 
fect the subsequent image subtraction. However, probably 
because our reference epoch was reduced in such a man- 
ner, we found that the optimal image subtraction, with the 
least numbers of artificial residuals, was obtained by follow- 
ing the same procedure. The photometric calibration of the 
first-epoch images was done by K04, reaching a precision for 
the zero points of about 0.05 mag. We calibrated our images 
relative to the first epoch by comparing the photometry of 
all the objects detected in both epochs. The mean of the 
differences between the two measurements was taken to be 
the difference in zero points. The final images reach 3<r lim- 
its (defined as in K04) of R = 27.4, i = 27.1, and z = 26.3 
(5a limits of R = 26.9, i' = 26.7), and z = 26.0), shallower 
than the reference image by 0.5 mag (R), 0.4 mag (i'), and 
0.4 mag (z'). 

Reduction of the KPNO images was done using the 
IRAF package MSCRED, following the same steps as with 
the Subaru data, except that PSF degradation was not ap- 
plied. The final images reach 3cr limiting magnitudes corre- 
sponding approximately to R — 25 — 26 mag (5<r limits of 
R = 24.5 — 25.3 mag), depending on the epoch. The depth 
reached implies a S/N improvement in the filterless data, 
relative to _R-band imaging, of a factor near 2, for two of our 
epochs, and only a minor gain for the remaining two. HST 
images were reduced using the PyRAF script MultiDrizzle 
(Koekemoer et al. 2002). 

We performed PSF matching and image subtraction be- 
tween the new and reference Subaru images in all bands, 
and between the four KPNO epochs, using both the soft- 
ware ISIS (Alard & Lupton 1998; Alard 2000) and the Com- 
mon PSF Method (CPM; Gal- Yam et al. 2004a; Gal- Yam 
et al. 2007, in preparation). Briefly, ISIS minimizes a spa- 
tially variable convolution kernel that degrades the PSF of 
the better-seeing image to the PSF of the worse-seeing im- 
age. CPM convolves each image with the measured PSF of 
the other, and thus the final PSF of both images is worse 
than the initial PSF of either image, but kernel-matching, 
and its inherent uncertainties in the presence of noise and 
pixellation, are avoided. Both procedures produce output 
images with nominally identical PSFs, while introducing a 
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minimal amount of noise. While ISIS is robust and stable, 
we find that in some cases it tends to produce more subtrac- 
tion artefacts than CPM, mimicking SNe near the cores of 
bright galaxies. We have therefore used both codes. 

As a consequence of the dithering, the final images have 
a field of view of 0.31 deg 2 , but with a substantial region 
along the edges where the S/N is significantly smaller. This 
is due to the different effective exposure on the fringes of the 
field. We therefore trim down the difference image to 81% 
of its full size, and remain with a total differenced area of 
0.25 deg 2 (which is, in fact, the size of the Suprime-Cam 
field). 

2.2 Spectroscopy 

We obtained spectra of some of our SN host galaxies (see 
Ej4j> , and of several hundred random galaxies in the SDF. Ob- 
servations were carried out on 2007 January 12, January 22, 
March 16, and April 12, using the Low-Resolution Imag- 
ing Spectrometer (LRIS; Oke et al. 1995) on the Keck I 
telescope for the first night, and the Deep Imaging Multi- 
Object Spectrograph (DEIMOS; Faber et al. 2003) on the 
Keck II telescope for the three other nights. Eight unique 
object masks and positions where used, with up to 12 host 
galaxies per mask (including also targets from Subaru ob- 
servations in February 2007, which are not analysed in this 
paper), and several tens of random galaxies per mask. Each 
mask was typically observed for 3 x 1800 s. 

The LRIS mask was observed with the 600 line mm 
grism blazed at 4000 A, and the 400 line mm -1 grating 
blazed at 8500 A, together with the D560 dichroic. This 
typically yields a wavelength range of ~3000-10,000 A, but 
the specific position of the slit on the mask shifts the wave- 
length range for each individual spectrum. The spectra have 
resolutions of ~3.5 A and ~7.0 A for the blue and red sides, 
respectively. 

All but one of the DEIMOS masks were observed with 
the 600 line mm -1 grating blazed at 7500 A, and the GG495 
order-blocking filter. We chose a wavelength range of ~5000- 
10,000 A, with the precise limits depending on each individ- 
ual spectrum (different slit positions). One mask was de- 
signed to target 1.5 < z < 2 galaxies as well as slightly 
higher- redshift SN hosts (see i )4.2|) . For this mask, we ob- 
served with the 600 line mm -1 grating and the OG550 order- 
blocking filter, with a central wavelength of about 8000 A. 
This yielded a typical wavelength range of 5500-10,500 A, 
without second-order contamination at longer wavelengths, 
due to the different filter. 

The 600 line mm -1 grating yields a full width at 
half-maximum intensity (FWHM) resolution of ~1.5 A, or 
~150 km s _1 , at 7500 A. This resolution is sufficient to re- 
solve many night-sky lines and the [O n] AA3726, 3729 dou- 
blet. Resolving night-sky lines is useful for finding emission 
lines in the reddest part of the spectrum, where sky lines 
are blended in low-resolution spectra. Resolving the [O n] 
doublet allows us to confidently identify an object's redshift 
even with a "single" line. 

The DEIMOS data were reduced using a modified ver- 
sion of the DEEP2 data reduction pipelinq 2 ], which bias cor- 

2 http:/ /astro. berkeley.edu/~cooper/deep/spec2d/ 



rects, flattens, rectifies, and sky subtracts the data before ex- 
tracting a spectrum (Foley et al. 2007) . The LRIS data were 
reduced using a combination of typical IRAF techniques and 
our own IDL procedures (Foley et al. 2007). The wavelength 
solutions were derived by low-order polynomial fits to the 
lamp spectral lines, and shifted to match night-sky lines at 
the positions of the objects. Standard-star spectra were ob- 
tained through a long slit on the same night, and were used 
to flux calibrate the spectra and remove telluric absorption 
(Matheson et al. 2001). To obtain the proper absolute flux 
scale and to correct for minor continuum differences due to 
slit losses, the spectra were scaled to match the i?-band and 
z'-band photometry of the galaxies. 

3 SUPERNOVA CANDIDATES 

3.1 Candidate Selection 

The z'-band difference image obtained using ISIS was 
scanned by eye to search for variable and transient objects in 
the new image. Morokuma et al. (in preparation) will present 
a study of the various variable objects they detect within the 
subframes that constitute our reference image set. They have 
also identified nearly a thousand AGNs in the SDF based 
on their long-term i'-band variability. In our search, these 
galaxies were therefore ignored (clearly non-nuclear variable 
sources in these objects were searched for, but none was 
found). The remaining variable candidates were examined 
as follows, in order to reject other non-SNe. 

(i) We compared the z'-band difference images relative to 
the reference SDF epoch for the first and second nights of the 
data separately, to identify and reject moving objects (as- 
teroids) and possible subtraction artefacts and noise peaks. 

(ii) ISIS and CPM subtraction images were compared to 
test for possible subtraction artefacts, especially near bright 
galaxies. Candidates that had distorted or otherwise suspi- 
cious shapes in the ISIS subtraction image, and were clearly 
absent from the CPM image, were rejected. 

(iii) KPNO subtraction images were used to identify 
AGNs based on their non-SN-like (e.g., slowly rising) light 
curves. For 10 of our candidates, at the bright end of the 
distribution, we have detections also in one or more of the 
KPNO images, with a photometric behavior consistent with 
SNe, i.e., declining on time scales of a month or two. 

(iv) For every candidate found in the z' band, subtrac- 
tion images in the R and i' bands were also examined, and 
objects which showed suspect residual shapes, indicative of 
a subtraction artefact, were rejected. We note that no can- 
didate was rejected because of a non-detection in one band 
or another, since at least some high-z SNe are expected to 
be very faint or undetected in the observed-frame R and i' 
bands. 

We further note that the purity of our SN sample will in- 
crease retroactively as we obtain additional epochs of this 
field, and identify any remaining AGNs that may still con- 
taminate the sample. 

3.2 Detection Efficiency Simulations 

In order to assess our detection efficiency, artificial point 
sources, with characteristics matching as closely as possi- 
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Figure 1. Fraction of simulated SNe recovered as a function of 
z'-band magnitude. Error bars indicate 1<t Poisson uncertainties. 



ble those of the SN population, were planted blindly in the 
z'-band image and were searched for, together with the real 
transients. The simulated sample was constructed as follows. 
First, we calculated the photometric redshifts of the entire 
galaxy sample in the SDF (about 150,000 objects; see K04), 
using the Bayesian photometric redshift (photo-z) code BPZ 
(Bem'tez 2000). Briefly, BPZ fits a set of template galaxies 
to the photometry, with the possibility to use a prior on the 
redshift. We used the same 6 template galaxies as Bem'tez 
(2000), from Coleman, Wu & Weedman (1980) and Kinney 
et al. (1996). As a prior we took the redshift distribution 
of galaxies in the Hubble Deep Field, which has a depth 
comparable to that of the SDF (Fernandez-Soto, Lanzetta 
& Yahil 1999). The errors given by BPZ were used to se- 
lect a subsample with relatively robust redshifts (i.e., small 
nominal uncertainties), while requiring the selected galax- 
ies to have a redshift distribution similar to the one found 
in the Hubble Deep Field. A discussion of the photo-z de- 
termination for the actual SN host galaxies, along with our 
spectroscopic training set, are presented in Sj4] 

In order to account for the population of galaxies which 
are below our limiting magnitude, but can still host SNe, 
we extrapolated, in every redshift bin, the galaxy popula- 
tion using a Schechter (1976) luminosity function and as- 
signed those extrapolated galaxies random positions within 
our field. SN rates are expected to roughly follow the lu- 
minosities of the galaxies. We approximated the luminosity 
of each galaxy using its magnitude and luminosity distance 
(assuming here, and throughout the paper, cosmological pa- 
rameters H = 70 km s" 1 Mpc~\ Q. m = 0.3, SIa = 0.7). 
These luminosities were used to weight each galaxy's proba- 
bility to host a SN. The position of the simulated SNe within 
each host galaxy was again set to follow the light. We de- 
fined for each galaxy, using SExtractor (Bertin & Arnouts 
1996), 10 annuli, each containing one-tenth of the total flux 
of the galaxy. The radial position of the simulated SNe was 
randomly drawn to lie inside one of those 10 annuli. While 
the annuli are circular, the hosts are generally not. To ap- 
proximately compensate for this, the position angles of the 
SNe, relative to the galaxies' major axes, were chosen from a 
linear distribution, where the major axis position angle has 
the highest probability, and the orthogonal angle the lowest. 



Each fake SN was assigned a random magnitude, but 
was kept in the sample only if it was no brighter than the 
magnitude a SN la would reach at peak at the given redshift 
of the host. We thus generated a list of positions and magni- 
tudes for more than a thousand fake SNe. The objects were 
planted in the image using the IRAF task mkobject with 
a Gaussian profile. During the search, we could not distin- 
guish between real and fake SNe. Our resulting efficiency as 
a function of magnitude can be seen in Figure [T] Our re- 
covery rate is stable and nearly perfect, averaging 96%, up 
to z' = 25 mag, where it starts to decline. The efficiency 
is 50% at z' = 26.3 mag, and effectively reaches zero at 
z' = 27 mag. More than 80% of the 1000 simulated SNe 
are in the interesting magnitude range, between 25 and 27, 
where a good sampling of the efficiency curve is important. 
We have tested the effect of replacing the Gaussian profile 
of the fake SNe by an empirical profile measured from the 
image, using different template stars for different portions of 
the image, and found no changes in the efficiency compared 
to our previous simulations. 

3.3 Supernova Sample 

We found a total of 33 SNe, shown in Figure [2] with mag- 
nitudes in the range z = 23.7 to z' = 26.3. Table[T]lists the 
SNe and their properties. Our SNe do not have spectroscopic 
confirmation (indeed, spectroscopy for most of them is im- 
possible with existing telescopes) or light curves, therefore 
they do not satisfy the International Astronomical Union cri- 
teria for a "standard" appellation. Nonetheless, since we are 
confident in our search and screening process (and also for 
the sake of brevity), we will refer to them as SNe, rather than 
SN candidates. We denote the SNe from this March 2005 
run as "SNSDF0503.XX" , XX being a serial number, or- 
dered roughly according to the SN z'-band magnitude. The 
respective host galaxies are referred to as "hSDF0503.XX" . 
Ten of the 11 brightest SNe are also visible in the first im- 
age from KPNO. Apart from the 33 SNe, we detect several 
tens of candidates at fainter magnitudes, as we would expect 
based on our efficiency simulations (see 13. 2\ , but these are 
lcr to 2ct detections, which are mixed with an unknown num- 
ber of false positives that arise from random peaks in the 
noise distribution and from subtraction artefacts. We there- 
fore limit our sample to z' ^ 26.3 mag, which leaves only 
candidates that we judge to be unambiguously real. Coinci- 
dentally, z' ^ 26.3 mag is also our 50% detection efficiency 
limit. 

Aperture photometry of the SNe was performed on the 
difference R, i', and z' images, using SExtractor, and adopt- 
ing a 2" diameter circular aperture. To estimate the photo- 
metric uncertainty, we measured the magnitudes of the 1000 
simulated transients planted in the image for the efficiency 
measurements (see £|3.2[) . and took the root mean square 
(rms) in each magnitude bin to be the minimum photometric 
error for objects of that magnitude. The adopted uncertainty 
for each SN was taken to be the larger among the formal er- 
ror computed by SExtractor and the statistical error for the 
given magnitude from the simulated objects. For each SN 
we derived the local S/N by comparing the SN counts in 
the photometric aperture, to the standard deviation of the 
total counts in tens of identical apertures on nearby, sim- 
ilar, background. We also measure the offset between the 
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simulated SN positions as input to the images, and as found 
by SExtractor. All the simulated SNe are found within O'.'IO 
of their intended positions, and 93% are within 0'.'04. We 
use these results to estimate the accuracy to which an offset 
between a SN and its host-galaxy centre can be detected. 

The two SNe closest to their respective host centres 
have measured offset of 0'.'07 and 0'.'08, while all others have 
offsets greater than O'.'IO, effectively ruling out an AGN clas- 
sification for most, if not all, objects. 



4 SUPERNOVA HOST GALAXIES 

In this section, we determine which galaxy hosts every SN, 
and measure its properties. This allows us to measure the 
redshift of the SNe, and eventually, will permit a study of 
the correlations between the properties of SNe and their host 
galaxies. 



4.1 Identification and Photometry 

For our list of SNe, we compiled a list of potential host 
galaxies. We measured the Petrosian magnitudes (Petrosian 
f976) of the host galaxies on the reference (deeper) epoch, in 
seven bands, with SExtractor. Petrosian magnitudes allow 
one to measure the flux of resolved objects within a given 
fraction of the light profile, without the dependence on the 
amplitude of the surface brightness profile associated with 
"isophotal" magnitudes (e.g., Blanton et al. 2001; Graham 
et al. 2005). As for the SN photometry, we estimated the 
uncertainty in each magnitude bin using artificial sources 
with galactic profiles that we planted in the images. 

The hosts of the SNe were chosen based on the smallest 
separations between each SN and its neighboring galaxies, 
in terms of each of the galaxies' half-light radii, as measured 
with SExtractor in the i' band. For most of the SNe, the 
choice was obvious, since the SNe and their hosts were un- 
resolved from each other. 

For two SNe (SNSDF0503.15 and SNSDF0503.18) we 
found no plausible host within 6 half-light radii. Assuming 
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a Sersic (1968) model for the galaxy radial profile between 
n = 4 (de Vaucouleurs law; de Vaucouleurs 1948; Peng et al. 
2002) and n = 1 (exponential disk; Freeman 1970; Peng 
et al. 2002), we find that between 91% and 99.99% of the 
light, respectively, is within 6 half-light radii . The proba- 
bility that a SN belongs to a host at a distance beyond this 
limit is thus small. Due to the high density of sources in the 
SDF, using a larger limit would lead to confusion regarding 
the hosts of most of the SNe, including those SNe which are 
unresolved from a particular galaxy. The two SNe which we 
label as "hostless" are probably the result of a host that is 
below our detection threshold in all the photometric bands 
(though these SNe could, in fact, be hostless; see Gal- Yam 
et al. 2003). Alternative explanations are that these can- 
didates are actually flaring Galactic M-dwarfs or variable 
AGNs in faint hosts. The non-detection of a source at this 
position in the reference epoch, down to i' w 27.5 mag, 
means that a star with a quiescent absolute magnitude of 
Mi w 8 (e.g., The, Steenman & Alcaino 1984) would be at 
a distance > 80 kpc, making this option highly improbable. 
The possibility of a high-z AGN cannot be excluded. Nev- 
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ertheless, the x 2 criterion we apply to the candidates (see 
fj5| does not support this interpretation. Table [2] lists the 
SN hosts and their properties. 



4.2 Host Redshifts 

We find the redshifts of the host galaxies using a com- 
bination of photometry and spectroscopy. Out of 19 host 
galaxies with R < 25 mag, the approximate practical 
limit for spectroscopy, we observed 18, and managed to 
measure 16 redshifts. We also measured redshifts for a 
fainter host, hSDF0503.03, with R = 25.8 mag. The spec- 
tra are shown in Figure [3] The spectra of hSDF0503.16 
and hSDF0503.17, which were presumably bright enough for 
Keck (R — 23.56 mag and R = 24.91 mag, respectively), did 
not yield a conclusive redshift, probably due to slit imperfec- 
tions that impair proper sky-line subtraction. We note that 
no spectrum shows an AGN signature (except hSDF0503.35 
that hosts a LINER offset from the SN, see Sj5}. 

The highest spectroscopic redshift we obtained is z — 
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Figure 2. Image sections, 10" on a side, for the SNe in our sample. North is up, east is to the left. For each SN we show the reference 
and new epoch z'-band images (left and centre, respectively), and the difference image (right). The greyscale for all the difference images 
is the same, while for the other images it is chosen to emphasize host-galaxy details. 



1.564, for hSDF0503.40, slightly higher than the redshift 
of the host of SN 2003ak found by GOODS at z = 1.551 
(Riess et al. 2004a). This makes SNSDF0503.40 the SN with 
the highest host spectroscopic redshift reported to date. SN 
1997ff at 2 ~ 1.7, like other SNe in our sample having com- 
parable photometric redshifts, does not have a spectroscopic 
SN or host redshift (Riess et al. 2001). 

In order to assemble a training set for the derivation 
of photo-z, we measured 145 spectroscopic redshifts (123 of 
which are secure) for galaxies chosen randomly near the posi- 
tions of the SN hosts, at first with no selection criteria other 
than sufficient brightness for spectroscopy (i' < 25 mag) . In 
the last Keck run, we preferentially selected galaxies with 
photometric redshifts between 1.5 and 2, in order to bet- 
ter sample that region of parameter space. The redshifts for 
these 123 galaxies are listed in Table [3] in the online version 
of this paper. We supplement these data with spectroscopic 



redshifts for a sample of SDF galaxies that were obtained 
for other SDF-related projects (e.g., Kashikawa et al. 2003, 
2006; Shimasaku et al. 2006), usually selected for their excess 
flux in the narrow-band niters 7V.B816 and 7V.B921 (designed 
to detect Lyo? emission lines at redshifts near 5.7 and 6.6, 
respectively). From this sample we select 196 galaxies with 
secure z < 2 redshifts. This results in a set of 319 galax- 
ies with secure redshifts that allows us to tune and test the 
photo-z determination. While this sample is quite extensive, 
it is probably biased toward galaxies with bright emission 
lines, for which spectroscopic redshifts are easier to measure. 

For the derivation of the SN host photometric redshifts, 
we have used the recent version of the Zurich Extragalac- 
tic Bayesian Redshift Estimator (ZEBRA; Feldmann et al. 
2006). Briefly, ZEBRA works as follows. First, the photome- 
try catalog is corrected for systematic errors by finding resid- 
uals in a first fit to a set of basic galaxy spectral templates. 
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Table 1. Supernovae Discovered in the Subaru Deep Field. a Signal-to-noise ratio of the SN, in the z'-band image. b Probability that 
object is a SN la, based on SN-ABC (P07). c Approximate absolute magnitudes in rest-frame B band, unless rest-frame B does not fall 
within the observed R,i' ,z' bands, in which case we list rest-frame U or V absolute magnitudes, at high or low redshift, respectively. d 
Fainter than 3<r limiting magnitude in this band. 



Next, the templates are corrected using a training set with 
spectroscopic redshifts. Finally, the redshifts are derived us- 
ing a Bayesian methodology, where the redshift and best- 
fitting template for each galaxy are found iteratively, while 
using the resulting distributions (both in redshift and tem- 
plate space) as priors. Our initial templates are the same 
ones used by Bem'tez (2000). For the correction of the pho- 
tometry, we measured all the objects in the SDF, following 
the procedure described in £|3. 3 1 for the photometry of the 
SN host galaxies. The six basic templates that are used in 
BPZ (for galaxy types elliptical, Sbc, Scd, irregular, and 
two types of starburst) were interpolated to produce three 
intermediate galaxy types between every two consecutive 
templates. The resulting 21 templates were then corrected 
using a partial (early) version of the training set described 
above, of about 200 galaxies. The correction was done in 
three redshift bins (z < 0.5, 0.5 < z < 0.9, z > 0.9). This 
resulted in 84 templates (the 21 original interpolated ones 
plus 21 for each redshift bin). The resulting posterior red- 
shift distribution of all the galaxies in the SDF was used as 
a prior for the SN host redshifts. We set the limit on the 
absolute rest-frame B-band magnitude of the galaxy to lie, 
conservatively, between —26 and —10. 

Figure|4]shows the results of using ZEBRA on our train- 
ing set. First, one can see that the photo- z determination is 
usually in good agreement with the spectroscopic redshift. 
On the other hand, photo- z values above 1.8 are many times 
unreliable for this training-set sample. Restricting ourselves 



to photometric redshifts smaller than 1.8, where 93% of our 
sample is located (and probably all of our SNe), leaves us 
with 296 galaxies, with a dispersion of cr z /(l + z) = 0.08 
after rejecting five 4er outliers (less than 2% of this sam- 
ple). The scatter for our sample of 17 SN host galaxies is 
somewhat smaller, u z /(l + z) = 0.06, without a single out- 
lier rejected. Application of the Kolmogorov-Smirnov test 
to the two samples of residuals indicates that the difference 
is not significant. 

Improvement to the photo- z derivation could come from 
an extension of the imaging data to the near-infrared, as 
the rest-frame 4000-8000 A range is shifted, at z ~ 1.8, into 
the J, H, and K bands. With existing wide-field near-IR 
imagers (e.g., WFCAM on the UKIRT, Casali et al. 2001), 
and certainly future ones (e.g., NEWFIRM on the NOAO 
telescopes, Autry et al. 2003), deep imaging in those bands 
of well-studied fields, such as the SDF, would be extremely 
useful for projects such as ours. 

The output from ZEBRA that we use is the red- 
shift probability distribution function (z-pdf) for each host 
galaxy, which is obtained by marginalizing the full posterior 
distribution over all templates. The resulting z-pdfs of our 
host sample are shown in Figure [S] The z-pdf can have a 
complex structure — it is not necessarily a single peaked, 
well-behaved, function. For example, the distribution for 
hSDF0503.03 is clearly bimodal, while that of hSDF0503.12 
contains very little information, with all redshifts almost 
equally likely. Using the full z-pdfs in the subsequent anal- 
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hSDF0503.03, z=0.89 




hSDF0503.05, z=0:59 




5500 6000 6500 7000 7500 8000 8500 9000 9500 

Observed wavelength [Ang] 



ID 


a (J2000) 


5 (J2000) 


R mag 


Redshift 


Instrument* 1 


1 


13:23:51.556 


+27:11:51.53 


24.3 


1.488 


4 


2 


13:23:40.027 


+27:13:01.35 


23.7 


0.287 


4 


3 


13:24:37.571 


+27:13:24.44 


23.4 


0.736 


2 


4 


13:24:21.800 


+27:13:22.62 


23.1 


0.530 


2 


5 


13:24:30.644 


+27:13:37.35 


25.1 


1.466 


2 



Table 3. Keck Spectroscopy of SDF Galaxies. The full table is available in the electronic version of the paper. a Instrument and 
observation date: 1 = LRIS, 12 Jan. 2007; 2 = DEIMOS, 21 Jan. 2007; 3 = DEIMOS, 16 Mar. 2007; 4 = DEIMOS, 12 Apr. 2007. 
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ysis allows us to take this uncertainty into account. Out 
of the 31 host z-pdfs, 23 have widths w z /(l + z) < 0.08, 
where w z is the standard deviation of the best-fitting Gaus- 
sian to the z-pdf, while only three have w z /(l + z) > 0.2 
(namely hSDF0503.12, hSDF0503.19, and hSDF0503.38). 
For 16 of the 17 host galaxies for which we have a spec- 
tral redshift, it is very close (usually nearly identical) to the 
photo-z, with Az/(1 + z) ^ 0.1, while for the remaining 
galaxy (hSDF0503.26) it differs by only Az/(l + z) = 0.15. 
For these 17 galaxies, we use in the subsequent analysis a 
Gaussian z-pdf centered on the measured spectroscopic red- 
shift, with a width w z — 0.01. For the two hostless SNe 
(SNSDF0503.15 and SNSDF0503.18), we use a flat redshift 



prior between and 2. Table[5]lists the measured properties 
of the SN host galaxies, including the best-fitting redshifts 
and template types. 



The region of the SN host galaxy hSDF0503.31 was im- 
aged with HST/ ACS. As seen in Figure [5] this host is a 
face-on spiral galaxy, consistent with the photo-z best-fitting 
template, with the SN clearly offset from the centre by 0'.'35. 
While this z = 1.08 galaxy seems to be associated with the 
brighter spiral seen in the figure, the photo-z of the bright 
galaxy is much smaller, z ~ 0.5. 
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Figure 3. Keck DEIMOS spectra for 17 of the SN host galaxies. For most galaxies, the redshift is securely determined by several emission 
or absorption lines. In hSDF0503.03, hSDF0503.26, and hSDF0503.40, the single line is double peaked, making the identification as [On] 
AA3726, 3729 quite robust. Regardless, any other assignment predicts the presence of other lines which we do not detect. Gaps in the 
spectra correspond to regions with strong night-sky line subtraction residuals. 



5 SUPERNOVA CLASSIFICATION 

The next step in our analysis is the classification of the SNe 
into types, based solely on our single-epoch photometry and 
the redshift information for the hosts. P07 have recently de- 
veloped for this purpose an algorithm, the SN Automatic 
Bayesian Classifier (SN-ABC). Briefly, the SN-ABC is a 
template-fitting routine that compares the SN magnitudes 



to synthetic magnitudes derived from a library of template 
SN spectra of different types, ages, redshifts, and extinc- 
tions. The input to SN-ABC is the photometry of the SN 
and the z-pdf of its host, which is used as a redshift prior. 
The output is the probability that the SN is of Type la, 
Pi a , compared to the probability it is a core-collapse event, 
Pec = 1 — Pia, a posterior «-pdf, and a \ 2 value. The \ 2 cri- 
terion permits rejecting candidates with colours unlike those 
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0.5 1 1.5 2 

Spectroscopic redshift 



Figure 4. Comparison of the spectroscopic redshifts and the ZE- 
BRA photometric redshifts for the random sample of 319 galaxies 
(grey crosses), and for 17 SN host galaxies (empty diamonds). 
Error bars indicate the 68.3% confidence interval of the 2-pdf 
for each galaxy. The scatter of the random sample with photo- 
z smaller than 1.8 is <r z /(l + z) = 0.08 after excluding five 4<r 
outliers. For the SN host sample the scatter is <r z /(l + z) = 0.06. 

of SNe, usually peculiar SNe, AGNs, or other transients or 
variables, that were not rejected in previous steps based on 
their variability characteristics. As shown in P07, more than 
half of the AGNs can be rejected in this manner. A SN with 
value of Pia > 0.5 is considered to be of Type la, while ob- 
jects with Pi a < 0.5 are classified as CC-SNe. As shown in 
P07, Pia also serves as a confidence estimator. The closer it 
is to unity (zero), the more secure is the classification as a 
SN la (CC-SN). 

P07 tested the algorithm on different real and simu- 
lated datasets. The SN-ABC successfully classifies 97% of 
the SNe la from the Supernova Legacy Survey (SNLS; Astier 
et al. 2006), and 85% of the Type II-P SNe. Similar success 
fractions were obtained for the GOODS SNe, and for a sim- 
ulated sample. These simulations also show that for Type 
Ib/c SNe, success fractions are smaller (about 75%), while 
Type Iln SNe are often misclassified as SNe la. The resulting 
purity of a SN sample depends on the details of the survey. 
Consequently, in !}S]we perform simulations with parameters 
tailored to our survey, and use them to "debias" our sample, 
and to determine the range of type distributions consistent 
with the data. 

Table Q] lists the SNe we have found, their redshifts, 
and their ABC classifications. Among the 33 SNe, 22 are 
classified as SNe la by the algorithm. Among these 22 SN la 
classifications, 16 have Pi a values higher than 0.7, while 6 of 
the 11 CC-SNe are classified with similar confidence, Pec > 
0.7. 

Two SNe, SNSDF0503.24 and SNSDF0503.35, have ex- 
treme x 2 values of 22.2 and 30.1, respectively. As discussed 
above, this could be indicative of a non-SN object, since, 
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Redshift 

Figure 5. Redshift probability density functions (z-pdfs) of the 
SN host galaxies, as calculated with ZEBRA, and used as pri- 
ors for the SN classification. Grey crosses mark the spectroscopic 
redshifts, when available. Notice the complex shapes of the z- 
pdfs, and the agreement in most cases between the most proba- 
ble photo- z and the spectroscopic redshift. All z-pdfs have a total 
probability of 1; the vertical scale varies from frame to frame. 



as shown in P07, the SN-ABC assigns high v 2 values to 
most AGNs. However, both are well offset from the centres 
of their host galaxies. The host hSDF0503.35 is a "grand- 
design" face-on spiral galaxy (Fig. [2} having a LINER spec- 
trum (Ho, Filippenko & Sargent 1997) with an Ha luminos- 
ity of ~ 10 39 ergs -1 , typical of such objects. Thus, these are 
probably not AGNs (see, however, Gal- Yam 2005 and Gal- 
Yam et al. 2007, in preparation, who found a background 
AGN projected near a lower-redshift galaxy). Both of their 
respective hosts have relatively low (spectroscopic) redshifts, 
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1.5 < z < 2.0), using the peak of the z-pdfs as best-fitting 
redshifts. These "raw" distributions are shown in Fig.Qand 
listed in Table 0] 



Figure 6. Left: image section of HST-ACS image of the host 
galaxy of SNSDF0503.31. Centre: ij-band SuprimeCam image 
prior to SN explosion. Right: difference between SuprimeCam im- 
ages. There is a clear offset between the SN and its host centre, 
a face-on spiral galaxy. 



and are, in fact, the lowest-redshift SNe in our sample. Both 
are classified as CC-SNe, with high confidence, which in the 
context of the SN-ABC means that their colours resemble 
significantly more those of a SN II-P than a SN la. However, 
the large \ 2 values imply that neither SN type is a good fit. 
Since these are unlikely to be SNe la, we keep them in our 
CC-SN sample. Nevertheless, we will continue to monitor 
these positions in future observations, in order to test the 
background AGN option. 

We also note that hSDF0503.35 is classified by ZEBRA 
as an elliptical galaxy, in contradiction with its spiral mor- 
phology. This is probably due to the photometry being dom- 
inated by the bulge, while the SN could be in the disc. An- 
other CC-SN host, hSDF0503.23, is also classified as an ellip- 
tical, suggesting that, here as well, the classification, either 
of the host galaxy, or of its SN, is erroneous. However, we 
expect a fraction of the SNe to be misclassified, and we deal 
explicitly with the resulting uncertainties and biases in £0 

The posterior redshift, derived from the combined fit 
of each SN and its host, is generally quite close to the prior 
(host-galaxy-only based) redshift, but as shown in P07, these 
posterior redshifts can sometimes be biased, and therefore 
we generally do not use them in the next step. A clear ex- 
ample of disagreement between a prior and a posterior red- 
shift is SNSDF0503.38, which has a prior redshift of 1.87, 
while the SN photometry combined with this prior yields 
a classification as a CC-SN at z = 0.3. This apparent dis- 
crepancy is easily understood if one notices that the z-pdf 
for the host galaxy of this SN is very wide and unconstrain- 
ing, and hence the best-fitting redshift for the host is quite 
meaningless. For this object, as well as for SNSDF0503.12 
and SNSDF0503.19 that have similarly wide z-pdfs, and the 
hostless SNSDF0503.15 and SNSDF0503.18, we use the pos- 
terior redshifts that take into account the SN properties. 

For a rough idea of the properties of the SNe, we use 
the peaks of the z-pdfs as best-fitting redshifts, and find ap- 
proximate rest-frame absolute magnitudes of the SNe. First, 
we fit the de-redshifted fluxes with a linear function, to de- 
scribe the spectral energy distribution (SED) of the SN. We 
then perform synthetic photometry in the B band, unless 
rest-frame B does not fall within the observed R, i' , or z 
bands. In such a case, we derive rest-frame U or V abso- 
lute magnitudes, at high or low redshift, respectively. The 
resulting absolute magnitudes, listed in Table [TJ range be- 
tween — 13.2 and —17.3 for the CC-SNe, and between —15.1 
and -19.5 for the SNe la. 

We divide the samples of SNe la and CC-SNe into four 
redshift bins (z < 0.5, 0.5 < z < 1.0, 1.0 < z < 1.5, and 



6 UNCERTAINTIES, SAMPLE 

CONTAMINATION, AND DERIVATION OF 
INTRINSIC SUPERNOVA TYPE AND 
REDSHIFT DISTRIBUTIONS 

The ABC classification can introduce biases due to the de- 
pendence of the success fraction on the intrinsic parameters 
(type, age, redshift) of the SNe. We therefore use Monte 
Carlo simulations to determine the allowed range, and the 
most likely "true" redshift and type distribution, of the SNe 
(as we would have observed if we had perfect spectroscopic 
data for all of the SNe and their hosts). 

Following P07, we have randomly generated 10, 000 SN 
spectra of each of the main SN types: la, Ibc, II-P, and Iln 
(see, e.g., Filippenko 1997). We calculate the synthetic mag- 
nitudes of the fake SNe using the SEDs from the spectral 
templates of Nugent, Kim & Perlmutter (2002) for SNe la, 
Ib/c, and II-P. Since the Nugent et al. (2002) spectra of 
SNe Iln are theoretical blackbody SEDs, for this type only 
we use the templates from Poznanski et al. (2002). Abso- 
lute magnitudes and their dispersions are taken from Dahlen 
et al. (2004). The scatter in colour within each type is ap- 
plied to the CC-SNe by adding an intrinsic, normally dis- 
tributed, noise with a standard deviation of a = 0.2 mag, the 
value used by Sullivan et al. (2006a). The SNe la we sim- 
ulate are assigned different stretch values s, following the 
method described by Sullivan et al. (2006a). We simulate 
a Gaussian distribution of stretch values with an average 
of s = 1 and a dispersion of a = 0.25, truncated outside 
the range 0.6 ^ s ?C 1.4. We model the stretch-luminosity 
relation using the formalism Mbc = Mb — q(s — 1) (Perlmut- 
ter et al. 1999), where Mbc and Mb are the corrected and 
uncorrected B-band absolute magnitudes, respectively, and 
the correlation factor is a — 1.47. We also apply colour- 
stretch corrections using the method presented by Knop 
et al. (2003), by dividing the template spectra of normal, 
s = 1, SNe la, by smooth spline functions, in order to match 
their rest-frame UBVRI colours to those of SNe with vari- 
ous stretch values. 

To every simulated SN we assign an epoch, a host- 
galaxy extinction, and a redshift. The epoch is drawn from a 
uniform distribution, while the host extinction, Av, is drawn 
from the positive side of a Gaussian distribution with max- 
imum at zero, a = 0.2 mag for SNe la, and a = 0.5 mag 
for CC-SNe, truncated at Ay = 1 mag. The redshifts are 
drawn from the general galaxy population in the SDF us- 
ing the following scheme. A realistic photometric redshift 
is modeled, by drawing, for half of the SNe, the z-pdf of a 
random SDF galaxy with z < 2. To the other half of the 
simulated SNe we assign a spectroscopic redshift, by using a 
a — 0.01 Gaussian z-pdf. This mimics the redshift determi- 
nation characteristics of the real sample. We have measured 
the widths of the simulated sample's z-pdfs, as well as of the 
real host-galaxy sample, by fitting them with Gaussians, and 
find that the simulated and real samples have similar dis- 
tributions of widths. We draw a specific redshift for each 
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simulated SN from its z-pdf, and use it to redshift the sim- 
ulated spectrurrjf]. 

We then fold the simulated SN spectra through the ob- 
servational setup of our SDF observations, with their three 
bandpasses, photometric errors, and limiting magnitudes. 
We have measured the mean photometric errors of the SDF 
sample in each band as a function of magnitude, and assum- 
ing that the noise is normally distributed, have added it to 
each object. As with the real sample, we keep only objects 
with z'-band magnitudes brighter than 26.3, leaving 4, 000 
of each type. The simulated SNe are then blindly classified 
by the ABC into SNe la or CC-SNe, producing a mapping 
of success fraction vs. redshift for each SN type. We find the 
success fractions of the ABC in four redshift bins [z < 0.5, 
0.5 < z < 1.0, 1.0 < z < 1.5, and 1.5 < z < 2.0) for each of 
the four SN types. 

The success fractions are applied to the set of all possi- 
ble intrinsic type distributions (e.g., 30% SN la, 30% SN II- 
P, 20% SN Ib/c, and 10% SN Iln). Using steps of 2.5%, there 
are 12, 341 possibilities to distribute 100% of the observed 
SNe, in each redshift bin, among the four types. We com- 
pute, for each of these possibilities, the resulting fraction 
of objects that are labelled as SNe la in the sample (i.e., 
the sum of the fraction of SNe la that were correctly classi- 
fied), and of the fraction of CC-SNe that were misclassified 
as SNe la, as obtained by the average success fractions for 
each type of SN, in each redshift bin. Next, for each scenario, 
in each redshift bin, we find the binomial probability to draw 
the observed raw number of SNe la out of the total number 
of SNe in the bin, given the expected SN la fraction for that 
scenario. Since, in this scheme, the SN la fractions do not 
have equal prior probability (e.g., there is only one model 
with 100% SNe la, while there are many with 50% SNe la 
and various fractions for the CC-SN subtypes), we multiply 
each probability by a weight function that is inversely pro- 
portional to the number of scenarios with the same SN la 
fraction. Finally, we marginalize over all the scenarios in 
each redshift bin, and obtain the probability density func- 
tion for the true SN la fraction in that bin. From the prob- 
ability distribution function of the SN la fraction in each 
bin, we extract the most probable value, at the peak of the 
distribution, and the la uncertainties defined as the region 
that includes 68.3% of the probability density. We add to 
this error the corresponding fraction of the la Poisson un- 
certainty in that bin — that is, the Poisson uncertainty for 
the total number of SNe in the bin multiplied by the relevant 
fraction of SNe la. The most probable SN la fraction comes 
out the same as the SN la fraction given by the SN-ABC in 
two out of the four redshift bins. Only in the second redshift 
bin, 0.5 < 2 < 1.0, is the debiased most-probable fraction 
of SNe la reduced from 0.64 to 0.39. The raw and debiased 
redshift distributions of the SNe la and CC-SNe are listed 
in Table 2J and shown in Figure [7] We discuss the results in 
the following section. 



3 This is the most straightforward way to model reality, where 
the order is reversed, i.e., an object is at a particular redshift, and 
one measures its z-pdf. 
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Figure 7. Observed, raw (circles; slightly offset to the right for 
clarity), and debiased (lines), redshift distributions of SDF SNe la 
and CC-SNe, compared to the GOODS sample (squares; offset 
to the left). SDF error bars mark lcr classification and Poisson 
errors, added in quadrature. GOODS error bars denote la Poisson 
uncertainties. The debiased SDF and GOODS distributions are 
consistent to within the uncertainties, but there is a suggestion 
of more high-z SNe la in the SDF. 



7 COMPARISON TO THE GOODS 
SUPERNOVA SAMPLE 

Before calculating rates, which requires some assumptions 
and modeling of the SN properties, we compare our SN sam- 
ple to the most similar one in terms of depth and numbers, 
the GOODS SN sample. As part of the GOODS project, two 
fields of about 150 arcmin 2 each were observed with the ACS 
camera on HST every 45 days, over five epochs, in order to 
search for SNe. We have compiled from Riess et al. (2004a) 
and Strolger et al. (2004) all available photometry for the 42 
SNe found in the GOODS fields. These SNe cover a redshift 
range of 0.2 to 1.55, with a median of 0.76. The SNe have 
partial light-curve coverage, and not all have spectroscopic 
redshifts or types. However, GOODS is a complete sample, 
without any subsequent selection. 

Our survey and GOODS have different effective areas 
and observing strategies, for which we need to account prior 
to comparison. In terms of bands, the surveys are similar. 
The filter-plus-system bandpasses of the Suprime-Cam z fil- 
ter and of the FSTF850LP filter (used by GOODS) have ef- 
fective wavelengths at 9067 A and 9082 A, respectively. The 
GOODS detection efficiency function, as described by Strol- 
ger et al. (2004), is slightly deeper than our own, by about 
0.1-0.2 mag. We consider only SNe brighter than our 50% 
efficiency limit, z = 26.3 mag. For z'-band magnitudes 
brighter than 26.3, both surveys have similar efficiencies, 
and no correction is needed for a comparison to be made, 
if we exclude the fainter GOODS SNe. The SDF covers 
~ 900 arcmin 2 , and has a single epoch, since we searched 
for SNe only in the March 2005 images. GOODS covered 
300 arcmin 2 , imaged at five epochs, each separated by about 
45 days, i.e., 2-4 weeks in the typical SN rest frames. This 
means that the epochs are not independent and cannot be 
simply summed, since most SNe are discovered early and are 



Supernovae in the SDF: An Initial Sample, and Type la Rate, Out to z rs 1.6 17 



Sub-sample z < 0.5 0.5 < z < 1.0 1.0 < Z < 1.5 1.5 < Z < 2.0 



SN la - 


(raw) 







9 


10 


3 


SN la - 


de-biased" 


o.o: 


-1.7+0.0 
-0.0-0.0 


c J.+4.2+1.9 
°'°-3.8-1.4 


10 o+ + 4 - 3 

lu - u -4. 4-3.1 


, n +0.0+2.9 
J - u -1.5-1.6 


CC-SN 


- (raw) 


6 




5 








CC-SN 


- de-biased a 


6.o: 


-0.0+3.6 
-1.7-2.4 


o tr+3.8+2.9 
°-°-4.2-2.3 


n n +4.4+0.0 
u - u -o. 0-0.0 


n n +i. 5+0.0 
u - u -o. 0-0.0 


Total 




6 




14 


10 


3 



SN la rate 6 [10 -5 Mpc -3 yr _1 ] 0.0l^ 10.5^^ 8-lle.O 



Table 4. Supernova Redshift Distribution. "Errors are lcr classification and Poisson uncertainties, respectively. b SN la rate errors are 
classification and Poisson uncertainties added in quadrature 



detected in more than one epoch. From the photometry in 
Riess et al. (2004a) and Strolger et al. (2004), we find that, 
out of 42 SNe, 38 are detected (i.e., brighter than 26.3 mag, 
our SDF cutoff magnitude) in 1.45 epochs, on average, while 
the remaining 4 are too faint. Thus, the effective GOODS 
area is roughly 5 x 300 arcmin 2 /1.45 = 1036 arcmin 2 , about 
15% more than the SDF. Hence, for a direct comparison of 
the two samples, the number of GOODS SNe needs to be 
trimmed down from 38 to 38/1.15 = 33. Due to this reduc- 
tion factor, we will present non-integer numbers of SNe from 
GOODS in the discussion below. 

Thus, after the normalization for the different effective 
areas, there are (coincidentally) exactly the same number of 
SNe in our sample and in the "SDF-normalized" GOODS 
sample (33), 18.9 of which are SNe la (compared to the most 
probable value of 18.5 in the SDF debiased sample) and 14.1 
CC-SNe (most probable value 14.5 in the SDF). In Figure 
\7\ we compare the redshift distributions for the SN la and 
CC-SN subsamples from the two surveys. Since one of the 
CC-SNe in the GOODS sample (SN 2002fv) has no redshift 
information, we exclude it from this comparison. As can be 
seen, the two samples agree well in total size, in total type 
fractions, and in the CC-SN redshift distribution. 

However, the redshift distributions may differ for the 
SNe la. Our lowest-redshift SN la is at z « 0.5, while 
GOODS found SNe la down to z w 0.2. Furthermore, the 
SNe la in our sample reach higher redshifts, up to z ~ 1.6, 
rather than z ~ 1.4, in the "SDF-normalized" GOODS sam- 
ple. While formally the distributions are consistent, our data 
could indicate that the results in Strolger et al. (2004) and 
Dahlen et al. (2004) concerning the paucity of SNe la at high 
z, and the implications for the SN la delay time, may change 
once larger samples are obtained. On the other hand, one of 
our high-z SNe la, SNSDF0503.18, at z = 1.60, is hostless, 
and its redshift is based solely on the posterior fit from the 
SN-ABC to the SN photometry. If this redshift is in error, 
the number of SNe la in the highest-redshift bin would be 
more similar to that of GOODS. At z = 0.55, where the frac- 
tion of SNe la in our sample was most significantly reduced 
by the debiasing of our photometrically classified sample, 
lies the largest disagreement between competing measures 
of the SN la rates. Barris & Tonry (2006), who used pho- 
tometric methods to classify SNe, found a rate 4-5 times 
higher than measured by Pain et al. (2002) and Neill et al. 
(2006). Neill et al. (2006) have suggested that the sample 
of Barris & Tonry (2006) may have suffered from contami- 
nation by CC-SNe. Our raw sample of SNe la in that same 



bin, without debiasing, would certainly be contaminated by 
CC-SNe. 



8 THE TYPE IA SUPERNOVA RATE 

The SN redshift distributions that we have derived can be 
used to set constraints on the SFH and SN la delay func- 
tions. However, our SN sample is still small, and considering 
the similarity in the redshift distributions of the SDF and 
GOODS CC-SN samples, we postpone a CC-SN rate deriva- 
tion to future papers. Nevertheless, we calculate SN la rates, 
to examine how the differences in the distributions noted in 
SJ7] propagate into the underlying rates. We calculate the 
volumetric rate of SNe la, in four bins, 

_ Nla,j 

ria ^~ fr,(z)dV' 

where Ni a ,i is the number of SNe la in bin i, and rj(z) and 
dV are, respectively, the survey control time (see below) and 
the co- moving volume element (in Mpc~ 3 ), as a function of 
redshift z. The integration is done within the limits of each 
bin. The control time, rj(z), is found as follows. For every 
redshift z, we calculate the model z'-band SN la light curve, 
m z (t), and sum over the survey efficiency, e(m), at these 
magnitudes: 

v (z) = f«n.(t))%*. 

The effective redshift of each bin, < z > 4 , is defined as 
_ Jr)(z)zdV 

<z> >- fn(z)dv ' 

i.e., a mean redshift weighted by the volume and the control 
time. 

The uncertainties are the Poisson and classification un- 
certainties of the sample, added in quadrature. To investi- 
gate how errors in the efficiency function plotted in Fig. [T] 
might propagate to our rate measurement, we have fitted 
the efficiency curve with a Fermi-Dirac-like step-function, 
derived the uncertainties in its parameters, and have re- 
peated the rate calculation hundreds of times with randomly 
drawn efficiency curves. The resulting la dispersion in the 
SN rates is of the order of 10 _6 Mpc _3 yr _1 , i.e., about 30 
times smaller than the main sources of uncertainty. If, for 
some unknown reason, our detection efficiency curve is offset 
by 0.1 mag, the derived rates would change by about 15%, 
again less than the main sources of uncertainty. A signifi- 
cantly larger shift of the efficiency function,e.g., by 0.3 mag 
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toward bright magnitudes, would lower the efficiency to near 
zero at z' = 26.3 mag, and would be inconsistent with our 
having discovered SNe in this magnitude range (unless, of 
course, the intrinsic number of SNe rises sharply precisely 
at this brightness). We further neglect dust extinction in 
this rough derivation of the rate. As found by Riello & 
Patat (2005) and Neill et al. (2006), such an omission causes 
the rates to be underestimated by ~25% at most, which is 
smaller than the uncertainties already taken into account. 
We also do not attempt at this point to incorporate SNe 
with different stretch values into the control-time calcula- 
tion, which should have a negligible (if any) effect (e.g., Neill 
et al. 2006). The resulting rates are listed in Table [4] Fig- 
ure [8] shows the rates compared to other published results 
(as compiled by Neill et al. 2006) , with the SFR from Hop- 
kins & Beacom (2006) plotted to guide the eye. The SFR is 
scaled to fit the low-redshift (z < 0.5) measurements, where 
there are fewer conflicting estimations. 

In our lowest-redshift bin we have found no SNe la. Con- 
sequently, the error bar marks the upper limit arising from 
classification uncertainty only. Assuming the rate is in the 
range 10~ 4 ' 8 — 10~ 4 ' 2 yr _1 Mpc -3 , we expect between 1 and 
4 SNe la, given our survey volume and control time, which 
is consistent with us finding no SNe in this small-volume 
bin. In the 0.5 < z < 1 bin, our measurement appears to 
be inconsistent with the four highest-rate points of Barris & 
Tonry (2006), and of Dahlen et al. (2004). However, we note 
that these rate measurements show a conspicuous enhance- 
ment, relative to most of the other results, at 0.4 < z < 0.6. 
We further note that the error bars on the Barris & Tonry 
(2006) dataset are Poisson errors alone, and do not include 
any systematic uncertainties, such as classification errors 
that the sample may suffer from, at least to some extent. 
In the high-redshift bins, our rates are in good agreement 
with the Dahlen et al. (2004) results. Nevertheless, we find 
a ~50% higher rate of SNe la in the highest redshift bin. 
Overall, the central values of our measurements suggest a 
SN la rate that is fairly constant with redshift at high red- 
shift, and could be tracking the SFR. Such a behavior is 
expected at high redshift if the rate then is dominated by 
a prompt SN la channel with a short delay time, and con- 
trasts with the results of Dahlen et al. (2004), whose data 
suggest a declining SN la rate beyond z = 0.8 (Scannapieco 
& Bildsten 2005). 



9 CONCLUSIONS 

Large samples of high-z SNe are required to resolve crit- 
ical incompatibilities between recent results regarding the 
rates of SNe la, and their delay time from star formation 
to explosion. Here we have presented an initial sample of 
33 SNe from a new survey based on deep, single-epoch, re- 
imaging of the 0.25 deg 2 SDF. We have demonstrated that, 
for SN-rate purposes, such surveys can be carried out using 
ground-based observations, at a fraction of the cost of space- 
based data, and with the potential for large samples. Using 
Keck spectroscopy for more than half of the host galaxies, 
we have shown the reliability of photometric redshifts for 
such data, incidentally finding the SN with the highest host 
spectroscopic redshift reported to date, at z = 1.564. 

The SNe were classified using a novel Bayesian photo- 
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Figure 8. Type la SN rate from the SDF, compared to results 
from the literature (compiled by Neill et al. 2006). Empty circles 
are for Cappellaro et al. (1999), Hardin et al. (2000), Pain et 
al. (2002), Madgwick et al. (2003), Tonry et al. (2003), Blanc 
et al. (2004), and Neill et al. (2006); triangles are for Barris U 
Tonry (2006); filled circles are for Dahlen et al. (2004); and filled 
squares are for our derived SDF rate. The line shows the SFH 
from Hopkins & Beacom (2006), scaled to fit the z < 0.5 points. 

metric algorithm, using solely the SN photometry in three 
bands, and the host-redshift information, either photomet- 
ric or spectroscopic. The photometric classification of the 
SNe introduces biases for which we correct using simula- 
tions. The debiasing procedure we apply to the raw sample 
reduces significantly the fraction of SNe la in the redshift 
bin 0.5 < z < 1, which otherwise would have been heav- 
ily contaminated by CC-SNe. As discussed by Neill et al. 
(2006), such contamination may have affected the sample 
acquired by Barris & Tonry (2006), and could explain the 
discrepancy between their measurement of the SN la rate 
and the rates measured by Pain et al. (2002) and Neill et al. 
(2006). 

Our resulting sample is comparable to the GOODS SN 
sample, and is in good agreement with it, in terms of the to- 
tal numbers of SNe la and CC-SNe, and in terms of the red- 
shift distribution of the CC-SN sample. However, our sample 
of SNe la shows less of a decline at high redshifts than the 
GOODS sample does. This trend remains when comparing 
the SN la rates we derive to those derived by GOODS. We 
find at redshift z ~ 0.7 a rate which is similar to the values 
determined by Pain et al. (2002), Tonry et al. (2003), and 
Neill et al. (2006) at z ~ 0.5, but barely consistent with the 
results of Barris & Tonry (2006) and Dahlen et al. (2004). 
At z ~ 1.6, our rate determination is consistent with the 
result of Dahlen et al. (2004), but about 50% higher. 

These results, if confirmed by larger samples, may chal- 
lenge current conclusions, based on the paucity of SNe la 
found by GOODS at high redshift, concerning the SN la 
rate and the delay time. Additional epochs on this field, 
already being obtained, will enlarge our SN sample to the 
hundreds, and will test the reality of the apparent decline in 
the SN la rate at z > 1. Finally, our approaches to SN search 
and photometric classification are likely to become impor- 
tant in the era of synoptic telescopes such as Pan-STARRS 
(Kaiser et al. 2005) and the Large Synoptic Survey Telescope 
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(Stubbs, Sweeney & Tyson 2004), where the sheer numbers 
of SNe will require new approaches to analysis. 



ACKNOWLEDGMENTS 

We thank Nobuo Arimoto for his contribution to this 
project. We are also grateful to R. Feldmann, C. M. Carrollo, 
and P. Oesch for kindly providing us with the ZEBRA code, 
and for assisting us in its use. We wish to thank A. G. Riess, 
E. O. Ofek, J. D. Neill, and M. Way for helpful discussions 
and comments. This work was based on data collected at the 
Subaru Telescope, which is operated by the National Astro- 
nomical Observatory of Japan, and on observations collected 
at the Kitt Peak National Observatory (KPNO), National 
Optical Astronomical Observatories (NOAO), which is op- 
erated by AURA, Inc., under cooperative agreement with 
the National Science Foundation. Additional data presented 
here were obtained at the W. M. Keck Observatory, which 
is operated as a scientific partnership among the California 
Institute of Technology, the University of California, and the 
National Aeronautics and Space Administration; the Obser- 
vatory was made possible by the generous financial support 
of the W. M. Keck Foundation. The authors wish to recog- 
nize and acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most 
fortunate to have the opportunity to conduct observations 
from this mountain. We are grateful for support by the Inter- 
national Institute for Experimental Astrophysics at Tel- Aviv 
University. A.V.F.'s SN group at U.C. Berkeley is supported 
by US National Science Foundation grant AST-0607485, as 
well as by NASA/ffiST grant GO-10493 from the Space Tele- 
scope Science Institute, which is operated by AURA, Inc., 
under NASA contract NAS 5-26555. This research was sup- 
ported in part by the National Science Foundation under 
Grant No. PHY-0551164. 



REFERENCES 

Alard, C. 2000, A&AS, 144, 363 

Alard, C, & Lupton, R. H. 1998, ApJ, 503, 325 

Astier, P. et al. 2006, A&A, 447, 31 

Autry, R. G. et al. 2003, in SPIE Proa, Vol. 4841, 525 

Barris, B. J., & Tonry, J. L. 2006, ApJ, 637, 427 

Belczynski, K., Bulik, T., & Ruiter, A. J. 2005, ApJ, 629, 

915 

Benftez, N. 2000, ApJ, 536, 571 

Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 

Blanton, M. R. ct al. 2001, AJ, 121, 2358 

Cappellaro, E., Evans, R., & Turatto, M. 1999, A&A, 351, 

459 

Cappellaro, E. et al. 2005, A&A, 430, 83 
Casali, M. et al. 2001, in The New Era of Wide Field As- 
tronomy, ed. R. Clowes, A. Adamson, & G. Bromage (San 
Francisco: ASP, Conf. Ser. Vol. 232), 357 
Coleman, G. D., Wu, C.-C, & Weedman, D. W. 1980, 
ApJS, 43, 393 

Dahlen, T. et al. 2004, ApJ, 613, 189 

de Vaucouleurs, G. 1948, Annales dAstrophysique, 11, 247 
Faber, S. M. et al. 2003, in SPIE Proa, Vol. 4841, , 1657 



Feldmann, R. et al. 2006, MNRAS, 372, 565 
Fernandez-Soto, A., Lanzetta, K. M., & Yahil, A. 1999, 
ApJ, 513, 34 

Filippenko, A. V. 1997, ARA&A, 35, 309 

Foley, R. J., Smith, N., Ganeshalingam, M., Li, W., 

Chornock, R., & Filippenko, A. V. 2007, ApJ, 657, L105 

Forster, F., Wolf, C, Podsiadlowski, P., & Han, Z. 2006, 

MNRAS, 368, 1893 

Freeman, K. C. 1970, ApJ, 160, 811 

Gal- Yam, A. 2005, The Astronomer's Telegram, 586, 1 

Gal- Yam, A., & Maoz, D. 2004, MNRAS, 347, 942 

Gal- Yam, A., Maoz, D., Guhathakurta, P., & Filippenko, 

A. V. 2003, AJ, 125, 1087 

Gal- Yam, A., Maoz, D., & Sharon, K. 2002, MNRAS, 332, 
37 

Gal- Yam, A., Poznanski, D., Maoz, D., Filippenko, A. V., 

& Foley, R. J. 2004b, PASP, 116, 597 

Gal- Yam, A. et al. 2004a, ApJ, 609, L59 

Gal- Yam, A. et al. 2007, ApJ, 656, 372 

Graham, A. W., Driver, S. P., Petrosian, V., Conselice, 

C. J., Bershady, M. A., Crawford, S. M., & Goto, T. 2005, 

AJ, 130, 1535 

Greggio, L. 2005, A&A, 441, 1055 

Hillebrandt, W., & Niemeyer, J. C. 2000, ARA&A, 38, 191 
Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1997, 
ApJS, 112, 315 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 
Kaiser, N. et al. 2005, BAAS, 37, 1409 
Kashikawa, N. et al. 2004, PASJ, 56, 1011 (K04) 
— . 2003, AJ, 125, 53 
— . 2006, ApJ, 637, 631 

Kinney, A. L., Calzetti, D, Bohlin, R. C, McQuade, K., 
Storchi-Bergmann, T., & Schmitt, H. R. 1996, ApJ, 467, 
38 

Knop, R. A. et al. 2003, ApJ, 598, 102 

Kobayashi, C, Tsujimoto, T., & Nomoto, K. 2000, ApJ, 

539, 26 

Kobayashi, C, Tsujimoto, T., Nomoto, K., Hachisu, I., & 
Kato, M. 1998, ApJ, 503, L155 

Koekemoer, A. M., Fruchter, A. S., Hook, R. N., & Hack, 

W. 2002, in The 2002 HST Calibration Workshop, ed. 

S. Arribas, A. Koekemoer, & B. Whitmore, 337 

Li, W., Wang, X., Van Dyk, S. D., Cuillandre, J.-C, Foley, 

R. J., & Filippenko, A. V. 2007, ApJ, 661, 1013 

Madau, P., Delia Valle, M., & Panagia, N. 1998, MNRAS, 

297, L17 

Mannucci, F., Delia Valle, M., & Panagia, N. 2006, MN- 
RAS, 370, 773 

Mannucci, F., Delia Valle, M., Panagia, N., Cappellaro, E., 
Cresci, G, Maiolino, R., Petrosian, A., & Turatto, M. 2005, 
A&A, 433, 807 

Maoz, D., & Gal- Yam, A. 2004, MNRAS, 347, 951 
Matheson, T., Filippenko, A. V., Li, W., Leonard, D. C, 
& Shields, J. C. 2001, AJ, 121, 1648 

Muller, G. P., Reed, R., Armandroff, T., Boroson, T. A., 
& Jacoby, G. H. 1998, in SPIE Proa, Vol. 3355, , 577 
Neill, J. D. et al. 2006, AJ, 132, 1126 

Nugent, P., Kim, A., & Perlmutter, S. 2002, PASP, 114, 
803 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 
Oke, J. B. et al. 1995, PASP, 107, 375 
Ouchi, M. et al. 2004, ApJ, 611, 660 



20 Poznanski et al. 



Pain, R. et al. 2002, ApJ, 577, 120 

Panagi a, N., Delia Valle, M., & Mannucci, F. 2007, sub- 
mitted, |astro-ph/070 3409 

Peng, C. Y., Ho, L. C, Impey, C. D, & Rix, H.-W. 2002, 
AJ, 124, 266 

Perlmutter, S. et al. 1999, ApJ, 517, 565 
Petrosian, V. 1976, ApJ, 209, LI 

Poznanski, D., Gal- Yam, A., Maoz, D., Filippenko, A. V., 
Leonard, D. C, & Matheson, T. 2002, PASP, 114, 833 
Poznanski, D., Maoz, D., & Gal- Yam, A. 2007, AJ, ac- 
cepted, |astro-ph/0610129| (P07) 
Riello, M., & Patat, F. 2005, MNRAS, 362, 671 
Riess, A. G. et al. 2001, ApJ, 560, 49 
— . 2004a, ApJ, 607, 665 
— . 2004b, ApJ, 600, L163 

Scannapieco, E., & Bildsten, L. 2005, ApJ, 629, L85 
Schechter, P. 1976, ApJ, 203, 297 

Sersic, J. L. 1968, Atlas de Galaxias Australes (Cordoba, 

Argentina: Observatorio Astronomico, 1968) 

Sharon, K., Gal- Yam, A., Maoz, D., Filippenko, A. V., & 

Guhathakurta, P. 2007, ApJ, 660, 1165 

Shimasaku, K. et al. 2006, PASJ, 58, 313 

Strolger, L.-G. et al. 2004, ApJ, 613, 200 

Stubbs, C. W., Sweeney, D., & Tyson, J. A. 2004, BAAS, 

36, 1527 

Sullivan, M. et al. 2006a, AJ, 131, 960 
— . 2006b, ApJ, 648, 868 

The, P. S., Steenman, H. C, & Alcaino, G. 1984, A&A, 
132, 385 

Tody, D. 1986, in Proceedings of the SPIE, Vol. 627, , 733 
Tonry, J. L. et al. 2003, ApJ, 594, 1 

Yagi, M., Kashikawa, N., Sekiguchi, M., Doi, M., Yasuda, 
N., Shimasaku, K., & Okamura, S. 2002, AJ, 123, 66 
Yepes, G., Kates, R., Khokhlov, A., & Klypin, A. 1997, 
MNRAS, 284, 235 

Yungelson, L. R., & Livio, M. 1998, ApJ, 497, 168 



